New Variables for Classical and Quantum Gravity 
in all Dimensions III. Quantum Theory 

N. Bodendorfer^'2* T. Thietnann^'^J A. Thurn^^ 

^ Inst, for Theoretical Physics III, FAU Eriangen - NUrnberg, 
Staudtstr. 7, 91058 Eriangen, Germany 

^ Institute for Gravitation and the Cosmos Physics Department, 
Penn State, University Park, PA 16802, U.S.A. 

^ Perimeter Institute for Theoretical Physics, 
31 Caroline Street N, Waterloo, ON N2L 2Y5, Canada 

February 13, 2013 



Abstract 
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in our companion papers by Loop Quantum Gravity (LQG) methods. It turns out that all the 
tools prepared for LQG straightforwardly generalise to the new connection formulation in higher 
dimensions. The only new challenge is the simplicity constraint. While its "diagonal" components 
acting at edges of spin network functions are easily solved, its "off-diagonal" components acting 
at vertices are non trivial and require a more elaborate treatment. 
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1 Introduction 



In our companion papers [H [2] we developed the classical framework for a new connection 
formulation of General Relativity that is applicable in all spacetime dimensions D + 1 > 3. In 
3 + 1 dimensions, the current connection formulation is based on a triad and its corresponding 
spin connection. The miracle that happens in three spatial dimensions is that the defining 
representation of S0(3) is equivalent to its adjoint representation. Therefore, a connection 
and a triad carry the same number of degrees of freedom and can serve as a canonical pair 
on an extended phase space whose reduction by the SO (3) Gaufi constraint leads back to the 
ADM phase space. In order that the connection is Poisson commuting, a further miracle has 
to happen, namely the spin connection is integrable, i.e. can be obtained from a functional by 
functional derivation. These two miracles are reserved for D = 3. The observation that enables 
a connection formulation in higher dimensions as well is that the mismatch between the number 
of degrees of freedom of the D-bein and its spin connection can be accounted for by a new 
constraint in addition to the GauB constraint, which requires that the momentum conjugate to 
the connection comes from a D-bein. The details are a bit more complicated, we have to use 
SO{D + 1) rather than SO(-D), the D-hein is a generalised D-bein and the spin connection is a 
generalised hybrid connection, but this is the rough idea. 

The final picture is therefore a S0(D + 1) gauge theory subject to S0(D+1) Gaufi constraint, 
simplicity constraint, spatial diffeomorphism constraint and Hamiltonian constraint. Apart from 
the different gauge group which however is compact and the additional simplicity constraint, the 
situation is precisely the same as for LQG and the quantisation of our connection formulation is 
therefore in complete analogy with LQG. We can therefore simply follow any standard text on 
LQG such as [S*,"!] and follow all the quantisation steps. This way we arrive at the holonomy-flux 
algebra, its unique spatially diffeomorphism invariant state whose GNS data are the analogue 
for SO{D + 1) of the Ashtekar-Isham-Lewandowski Hilbert space, the analogue of spin network 
functions, kinematical geometrical operators such as the volume operator which is pivotal for 
the quantisation of the Hamiltonian constraint, the SO{D + 1) Gaufi constraint, the spatial 
diffeomorphism constraint, the Hamiltonian constraint and a corresponding Master constraint. 

The only structurally new ingredient is the simplicity constraint which constrains the type 
of allowed SO{D + 1) representations. When it acts at the interior point of edges, it requires 
that the corresponding SO{D + 1) representation is simple. However, when it acts at a vertex, 
the constraint splits into several linearly independent ones which are not mutually commuting 
and do not close on themselves. The situation here is similar to the situation in spin foam 
models O El [TJ El O [10] where similar constraints at the discretised level for S0(4) arise while 
ours are for SO{D + 1) in the continuum. We propose to solve these anomalous components of 
the simplicity constraints as in [H El [71 [8] by passing to a corresponding Master constraint and 
subtracting its spectral gapj^ 

The manuscript is organised as follows: 

In section two we define the SO{D + 1) holonomy-flux algebra and the corresponding Hilbert 
space representation. In section three we implement the kinematical constraints, that is Gaufi, 
simplicity and spatial diffeomorphism constraints. In section four we develop kinematical geo- 
metrical operators, specifically D-dimensional area and volume operators. Lower dimensional 
operators such as length operators etc. can be constructed similarly but are left for future publi- 
cation. Finally, in section five we quantise the Hamiltonian constraint. The presentation will be 

^The fact that the gauge group is compact makes sure that the spectrum of this Master constraint is pure 
point. 
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brief since all the constructions literally parallel those of LQG. We therefore refer the interested 
reader to [1] for all the missing details. 

2 Kinematical Hilbert Space 

The construction of the kinematical Hilbert has been performed in |1H [T2l [131 fl^ [T5l [TBJ for 
four and higher space-time dimension and arbitrary compact gauge group. These results apply 
for the case considered here, since we are using the compact group SO(D-|- 1) irrespective of the 
signature of the space-time metric. We therefore only cite the main results in this section and 
introduce notation needed later on. 

Since the Poisson brackets between Aaij and tt^^^ are singular, we have to smear them 
with test functions. In order to obtain non-distributional Poisson brackets, smearing has to 
be done at least D-dimensional in total. Aaij is a one-form, thus naturally smeared along a 
one-dimensional curve. From tt"'^'^ , being a vector density of weight one, we can construct the 
so{D -|- 1) - valued pseudo {D — l)-form {*TT)ai...ar,-i '■— t^""^^ (-aai...ao-i'''iJ which is integrated 
over a, {D — l)-dimensional surface in a background- independent way. These considerations lead 
to the definitions of holonomies and fluxes, which yield a natural starting point for a background 
independent quantisation. In the following, we choose (t/j)^l = \ {SfSjL — S^Sil) as a basis 
of the Lie algebra so(D -|- 1). 

2.1 Holonomies, Distributional Connections, Cylindrical Functions, Kine- 
matical Hilbert Space and Spin-Network States 

Denote by A the space of smooth connections over a. We define the holonomy hc{A) £ S0{D+1) 
of the connection A £ A along a curve c : [0, 1] — >• o" as the unique solution to the differential 
equation 

^^hcM) = hcM)A{c{s)), /ie„ = lD+i, hM) = hcAA), (2.1) 
where Cs{t) := c{st), s G [0, 1], A{c{s)) := A^^^ {c{s))tjjc°'{s). The solution is explicitly given by 

KiA) =Vexp( fA)= 1d+i + f][ dh [ dt2 . . . [ (it„^(c(ti)) . . . A{c{tn)), (2.2) 

where V denotes the path ordering symbol which orders the smallest path parameter to the left. 
Like in 3 -|- 1 dimensional LQG, we will restrict ourselves to piecewise analytic and compactly 
supported curves. 

The holonomies coordinatise the classical configuration space. In quantum field theory it 
is generic that the measure underlying the scalar product of the theory is supported on a 
distributional extension of the classical configuration space. For gravity, this enlargement of the 
configuration space is done by generalising the idea of a holonomy. Since the equations 

h,oc' {A) = hc{A)hc' {A) Vi (A) = hciA)-^ (2.3) 

hold, we see that an element ^ G ^ is a homomorphism from the set of piecewise analytic paths 
with compact support V into the gauge group. We now introduce the set A := IIom(P, SO{D + 1)) 
of all algebraic homomorphisms (without continuity assumptions) from V into the gauge group. 
This space A is called the space of distributional connections over a and constitutes the quantum 
configuration space. The algebra of cylindrical functions Cyl(^) on the space of distributional 
SO(D + 1) connections is chosen as the algebra of kinematical observables. The former algebra 



4 



can be written as the union of the set of functions of distributional connections defined on piece- 
wise analytic graphs 7, Cyl(^) = U^Cyl^{A) / ~. Cyl^(^) is defined as follows. A piecewise 
analytic graph 7 G u consists of analytic edg es ei,...,en, which meet at most at their endpoints, 
and vertices vi,...,Vm- We denote the edge and vertex set of 7 by E{'y) (|£'(7)| = n) and ^(7) 
(1^(7)1 = "^); respectively. A function f-y G Cyl^(^) is labelled by the graph 7 and typically 

looks like f^{A) = (^/lei(A), /ie|^| (A)) , where : SO{D + 1)1^1 C. One and the same 

cylindrical function / G Cyl(^) can be represented on different graphs leading to cylindrically 
equivalent representations of that function. It is understood in the above union that such func- 
tions are identified. We will denote the pullback of a function f-y defined on 7 on the biggeij^ 
graph 7' 7 via the cylindrical projections by Py^- Then, the equivalence relation just men- 
tioned can be made more explicit, ~ /y iff p*„.yf-y = PyiyfY V7, 7' -< 7". The pullback on 
the projective limit function space will be denoted by p*. The functions cylindrical with respect 
to a graph that are N times differentiable with respect to the standard differentiable structure 
on SO{D + 1) will be denoted by Cyl:^(:4) and Cyl^(:4) := U^Cyli^(:4)/ ~. 

Since in the end we are interested only in gauge invariant quantities, after solving the 
GauB constraint (classically oder quantum mechanically) we have to consider the algebra of 
cylindrical functions on the space of distributional connections modulo gauge transformations 
Cy\{A/Q). For representatives of elements / of this space, the complex- valued function F.y 
on SO{D + 1)1^1 has to be such that f^yiA) is gauge invariant. We will slightly abuse notation 
and use the same notation for the new projectors py-y : Ay/Gy — )• A^/Q-y. There is a unique 
\n\ [18] choice of a diffeomorphism invariant, faithful measure /io on A/G which equips us with 
a kinematical, gauge invariant Hilbert space := L2 (^A/Q, dfiQ^ appropriate for a representa- 
tion in which A is diagonal. This measure is entirely characterised by its cylindrical projections 
defined by 



[di^o{A)f{A) = [_diio,^{A)fyiA) 
JA/g JA/g 



SO(D+1)I^(t'I 



eeE{-y) 



F^{hi,...,h\E\), (2.4) 



where hh is the Haar probability measure on SO(D + 1). 

An orthonormal basis on Ti^ is given by spin- network states |19[ [20l I21j . which are defined 
as follows. Given a graph 7, label its edges e G E{'y) with non-trivial irreducible representations 
vTAg of S0(D + 1), i.e. Ae is the highest weight vector associated with e, and its vertices v G V{'y) 
with intertwiners Cy, i.e. matrices which contract all the matrices T^Ke{he) for e incident at v 
in a gauge invariant way. A spin-network state is simply a C°° cylindrical function on A/Q 
constructed on the above defined so-called spin-net, T^^^[A] := tr ®^.^iT^ K^^{hei{A)) ■ (^^^J^iCj , 
where A = (Ag), c = (c^) have indices corresponding to the edges and vertices of 7 respectively. 



2.2 (Electric) Fluxes and Flux Vector Fields 

Since ir""^^ are Lie algebra-valued vector densities of weight one, (*vr)aj. := t^"'^'^ €-aai...aD-i'^iJ 
is a pseudo [D — l)-form and is naturally integrated over a (Z) — l)-dimensional face S. We 
therefore define the (electric) fiuxes 

^"(5) := j n,j(*7r)" = / nijT:^'Jeaa,...ao-.dx^' A ... A dx^°-\ (2.5) 
Js Js 

^The graph 7 can be enlarged by e.g. adding or subdividing edges. See e.g. [Ij for a precise definition of the 
partial order on tame subgroupoids defined by graphs. 
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where n = n^'^Tjj denotes a Lie algebra- valued scalar function of compact support. We again 
restrict to piecewise analytic surfaces S, to ensure finiteness of the number of isolated intersection 
points of S with a piecewise analytic path. In order to compute Poisson brackets, we have to 
suitably regularise the holonomies and fluxes to objects smeared in D spatial dimensions. A 
possible regularisation in any dimension is given in |^. Removal of the regulator leads to the 
following action of the Hamiltonian vector fields Yn{S) corresponding to 7rn{S) on adapted 
representatives f^yg 

= E <^.S)n'\e^S)R\j^^,. (2.6) 

eeE(7s) 

/^g is an adapted representative of the cylindrical function / G Cyt^i^A) in the sense that all 
intersection points of S and 7s are beginning points 6(e) of edges e G -£^(75) (this can always 
be achieved by suitably splitting and inverting edges). In the above equation, e(e, S") is a type- 
indicator function, which is +(— )1 if the beginning segment of the edge e lies above (below) the 
surface S and zero otherwise. R\j {L'fj) is the right (left) invariant vector field on the copy of 
SO{D + 1) labelled by e, 

The algebra of right (left) invariant vector fields is given by 

/ 1 1 

Ru^Rkl] = ^^e,e'{VJKRlL+VlLRjK-VlKRjL-VJLR'lK), 

'Rh^Li^] = 0, (2.8) 
and analogously for L^j j- We remark that, in order to calculate functional derivatives, we had 



to restrict / to ^ in the beginning. The end result (2.6), however, can be extended to all of 
A. Following the standard treatment, these vector fields are generalised from adapted to non- 
adapted graphs and shown to yield a cylindrically consistent family of vector fields, thus they 
define a vector field Yn[S) on A. The Yn{S) are called flux vector fields. 

On the Hilbert space defined in section |2.1[ the elements of the classical holonomy-fiux 
algebra become operators which act by 

y„(5)-V := ihKPYn{Sm, (2.9) 

where the right hand side is the action of the vector field Yn{S) on the cylindrical function tp. 
The appearance of (3 is due to the fact that we defined the fluxes using tt, whereas the momenta 
conjugate to the connection is given by ^'^^vr = -gTr. The momentum operators 1^(5), with dense 
domain Cyl^, can be shown to be essentially self-adjoint operators on T-l^ analogously to the 
(3 + l)-dimensional case [13] . 



3 Implementation and Solution of the Kinematical Constraints 
3.1 Gaufi Constraint 

Working with the gauge invariant Hilbert space from the beginning, the Gaufi constraint is 
already solved. Yet we want to summarise its implementation on the gauge variant Hilbert 
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space % = L2 d^'o) , since we want to compute quantum commutators of the constraint with 
the simphcity constraint in the next section. The implementation (as well as the solution) of 
the Gaufi constraint can be copied from the (3 + l)-dimensional case without modification. 

According to the RAQ programme, we choose the dense subspace $ = Cyl°°{A) in the 
Hilbert space. Then, we are looking for an algebraic distribution L £ ^' such that the following 
equation holds 



Lipl 



ee-B(7); v=h{e) £££(7); v=f{e) 



/7 = (3.1) 



for any v G ^(7), any graph 7 and € Cyl'^{A). The general solution for L is given by a 
linear combination of (i/',.), where £ T~L^ is gauge invariant. Thus, for an adapted graph 7' 
(all edges outgoing from the vertex v in question), gauge invariance amounts to vanishing sum 
of all right invariant vector fields at a vertex, 

E = 0- (3-2) 

eeE(7'); v=b{e) 

3.2 Simplicity Constraint 
3.2.1 From Classical to Quantum 

Classically, vanishing of the simplicity constraints S^{x) = \^ijklm'^^^'^ {^)'^^^^{^) 
points X S 0" is completely equivalent to the vanishing of 

Cm(5"', -S""") := ^hm^ ^f^o-i)^'{D-if iJKLM'^^'^iS'^)'^^^iSd) (3-3) 

for all points x G a and all surfaces S^,S'^ C a containing x and shrinking to x as e, e' tend 
to zero. More precisely, we use faces of the form : (-1/2,1/2)^-1 ^ ^. {^^i, ...,ud-i) ^ 
S^{ui, ...,ud-i) with semi-analytic but at least once differentiable functions ^^(mi, md-i) and 



5^(0, ...,0) = X, and define 5f (ui, -ud-i) := S^{eui, euD-i)- We find that (2.5) becomes 
(with the choice njj = ^I^Sj^) 

1 If 

71^-^(5^) = / dui...duD~ieaai...aD-iidS'"'^/dui){ui,...,UD-l) X ■■■ 

e^^~^> J{~e/2,e/2)D-l 

x{dS^^^-'/duD-i){ui,...,UD-i) Tr'''-\S^Ui,...,UD-l)) 
= n,(5)^"^-^(x) + 0(e) (3.4) 

with na{S) = eaai...ao-i(55"'"V5wi)(0, ...,0) X ... X {dS'"'i^-^/duD-i)iO,.-.,0), from which the 
claim follows. Now, similar to the treatment of the area operator in section [4?H we just plug in 
the known quantisation of the electric fluxes and hope to get a well-defined constraint operator 



we find for a representative f-y^^, of / G Cyl^ {A) on a graph jss' adapted to both and S'^ 



in the end. Using the regularised action of the flux vector fields on cylindrical functions (2.6) 

1 



'€ 

,,l™0 eiD-l)^{D-lflJKLM Yl Y 



g(D-l)g'(D-l)^Af('S'^''5"')7ss'[-^7ss']- (3.5) 
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The flux vector fields only act locally on the intersection points e n 5, e G E{'yss')- Therefore, 
in the second line we used that for small surfaces Sf, , the action of the constraint will be 
trivial expect for x (and of course only non-trivial if x is in the range of 755'), thus independent 
of e. In the limit e, e' — )• the expression in the last line of the above calculation clearly diverges 

except for C/ = 0, where the whole expression vanishes identically. Since the kernels of the 

constraint operators C and C coincide, we can work with the latter and propose the constraint 
(omitting the ~ again) 



C''{S,S',x),p;f, = p;,,,e^™ e{e,Sne{e',S'nRhR'KLP;ssnfi 

e,e'e{e"e£;(7_gs,),6{e")=x} 
— on* JJKLM ( Tjup Tjdown\ ( ryup' r>down'\ f (o c\ 



where R!fP := Eeei?(7ss'),b(e)=x,^(e,5{'))=i ^/J similar for \ In the following, wih 

drop the superscript ^ for the surfaces for simplicity. 

The proof that the family (7*^(5, 5', x) is consistent and defines a vector field C^^ {S, S' ,x) 
on A follows from the consistency of Yn{S). To see that the operator is essentially self-adjoint, 
let ^ be the finite-dimensional Hilbert subspace of li? given by the closed linear span of spin 
network functions over 7 where all edges are labelled with the same irreducible representations 
given by vr, = (B^y^^Ti^^- Given any surfaces S, S' we can restrict the sum over graphs to 
adapted ones since we have C ^^^g, for the choice tt^, = vTe with E{'yss') 9 e' C e G 

E{'y). Since C^'\S, S', x) preserves each T-L^ ^, its restriction is a symmetric operator on a finite- 
dimensional Hilbert space, therefore self-adjoint. To see that it is symmetric, note that the right 



hand side of the first line of (3.6) consists of right-invariant vector fields which commute. This is 
obvious for the summands with vector fields acting on distinct edges e ^ e' , and for e = e' note 
that [-Rfj, Rkl\ antisymmetric in (U) o {KL) and thus vanishes if contracted with ^^J^lm ^ 
Now it is straight forward to see that C'^^ {S, S' , x) itself is essentially self-adjoint. 

Note that we did not follow the standard route to quantise operators, which would be to 
adjust the density weight of the simplicity constraint to be +1 (in its current form it is +2) and 
quantise it using the methods in |22]- Rather, the quantisation displayed above parallels the 



quantisation of the (square of the) area operator in 3+1 dimensions and indeed we could have 

considered J d^^^u ^J\n^Ji^~S^\ for arbitrary surfaces S and would have arrived at the above 

expression in the limit that S shrinks to a point without having to take away the regulator e (the 
dependence on two rather than one surface can be achieved, to some extent, by an appeal to the 
polarisation identity). If we would have quantised it using the standard route then it would be 



necessary to have access to the volume operator. We will see in section 4.2 that for the derivation 
of the volume operator in certain dimensions in the form we propose, which is a generalisation 
of the 3 + 1 dimensional treatment, we need the above simplicity constraint operator to cancel 
some unwanted terms. Of course, there might be other proposals for volume operators which 
can be defined in any dimension without using the simplicity constraint. Still, the quantisation 
of the simplicity constraint presented here will (1) give contact to the simplicity constraints used 
in spin foam models and (2) enable us to solve the constraint in any dimension when acting on 
edges. Its action on the vertices, i.e. the requirements on the intertwiners, is more subtle and 
we propose to treat it using the Master constraint method. We will first present the action on 
edges and afterwards derive a suitable Master constraint. For following calculations, note that 
we always can adapt a graph to a finite number of surfaces. Furthermore, it is understood that 
all surfaces intersect 7' in one point only (we may always shrink the surfaces until this is true). 
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3.2.2 Edge Constraints and their Solution 

The action of the quantum simpUcity constraint at an interior point x of an analytic edge 
e = ei o (62)"^ for both surfaces S, S' not containing e (otherwise the action is trivial) is given 

by 



C''{S,S',x)p*f- 



7^7 



2e' 



IJKLM ^ 
IJKLM 
IJKLM pei 
^^IJKLM j^ei pei 



^IJ-^KLP-ya^i-yf-y^ 



(3.7) 



where the sign is + if the orientation of the two surface 5, S' with respect to e coincides and — 
otherwise. In the second and fourth step we used gauge invariance at the vertex v of an adapted 



graph, Yl 



eG-B(7); v=b{e) ^IJ 



f- 



'yss' 



0, and in the third step we used that [i?"^^ , i?"^^] = 0. This 

(3.8) 



leads to the requirement on the generators of SO(D + 1) for all edges 

T[IJTKL] = 0. 

The so-called simple representations of S0(L' + 1) satisfying this constraint were classified in 
Irreducible simple representations are given by homogeneous harmonic polynomials ^j^^^^ of 
degree A^, in any dimension labelled by one positive integer N. In this sense, there is a similarity 
between the simple representations of SO{D + 1) and the representations of S0(3) (which all 
can be thought of as being simple). In particular, for + 1 = 4 we obtain the well-known simple 
representations of S0(4) used in spin foams labelled by = j = J • 

The commutator with gauge transformations at an interior point x of an analytic edge 
e = ei o (62)"^ (ei, 62 outgoing at x) yields, analogously to the classical calculation, 

G^gg, [A],^ {S,S ,x)^^^^, 

±A^^(x)e^™j^(i?-^ + R%) , {RTj - Rfj) {K, - m,)] 

D-3 



i=l 
D-3 

i=l 



(3.9) 



Two constraints acting at the same interior point x of an edge e = ei o (62) ^ commute weakly. 
Using the gauge invariance of C/ if / is gauge invariant, we find 



C^^iS, S',x),C^{S",S"',x' 



N ( Qll Qlll 



IJKLM OPQRN 



^IJ^KL^ ^OP^QR 



+ 0{Cf-,) + 0{Gf^ 



^ei ,rot 



*5^^^. (eR^' ■ C^i'™* + [(7^1.™*, ei?^i] + eR^' ■ C'^i-™* ) /, 



rot,rot 



/7 



(3.10) 
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which can be seen by the fact that the simphcity on an edge is quadratic in the rotation generator 
R^^ on that edge, and we used the notation 



Z)-3 

A^'^m; e^s^^^^i-^-i^.'^'+i-^«-3i?^^i?^^ =: A . C^.™* (3.II) 

i=l 



for a simphcity with a infinitesimal rotation acting on the muhi-index M (cf. (3.9)). Here, we 
chose a graph 7 adapted to all four surfaces S, S', S" , S'" . Note that classically, the Poisson 
bracket of two simplicity constraints vanishes strongly, whereas in the quantum theory this is 
only true in a weak sense. Still, the simplicity constraints acting on an edge are thus non- 
anomalous and can be solved by labelling all edges by simple representations of SO(D + 1). 

3.2.3 Vertex Master Constraint 

When acting on a node then, like the off-diagonal constraints in spin foam models, the simplicity 
constraints will not (weakly) commute anymore. Therefore, we are not allowed to introduce these 
constraints strongly and have the options of either trying to implement them weakly [5J or using 
a Master constraint. We will follow the latter route and give a proposal of how to construct a 
Master constraint of the simplicity constraints at the nodes. To reduce complexity, we try to 
find a both necessary and sufficient set of simple "building blocks" of the simplicity constraint 



at the node and construct a Master constraint using these. Considering (3.6), an obviously 
sufficient set of building blocks at the vertex v is given by 

R[ijRi<L]fi = Ve, e' e {e" G E{j); v = b{e")}. (3.12) 

For necessity, we have to prove that we can choose surfaces in such a way that these building 
blocks follow. Note that it has already been shown in p3] that all right invariant vector fields 
R'^ for single edges e can be generated by the Y{S), but the construction involves commutators 
of the fluxes. Since the simplicity constraints acting on vertices are anomalous, we cannot use 
commutators in our argument. Instead, we will construct the right invariant vector fields by 
using linear combinations of fluxes only. To this end, we will prove the following lemma: 

Lemma. 

For each edge e G E{v) at the vertex v we can always choose two surfaces S, S, such that the 
orientations with respect to S, S of all edges hut e coincide. 

The intuitive idea of how to find these surfaces is to start with a surface containing the edge 
e while intersecting all other edges e' G E{v), e' 7^ e transversally, and then slightly distort this 
surface in the two directions "above" and "below" defined by the surface, such that the edge 
e in consideration is once above and once below the surface, while the orientations of all other 
edges with respect to the surfaces remain unchanged, in particular none of them lies inside the 
surfaces. When subtracting the flux vector flelds corresponding to the two distorted surfaces, 
all terms will cancel except the terms involving R'^. 

Proof. To prove the statement above, two cases have to be distinguished: (a) the case where 
no e' S E{v) is (a segment of) the analytic extension through v of the edge e and (b) the case 
where e has a partner e which is a analytic extension of e through v. 

Case (a): The construction of the surface S^^e with the following properties 

1. Se C Sv,e for some beginning segment Se of e, and the other edges e' G E{v)., e' ^ e intersect 
Sv^e transversally in v. 
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2. For e' G E{v),e' / e: e' n 5^,e = -y, and for e' ^ ^(w), e' n S'j,,e = 0. 



is given in [21] and we summarise the result shortly. An analytic surface (edge) is completely 
determined by its germ [S]^ {[e]v) 

^ mi!...mz)_i! 

mi,...,mi3„i=U 

^(*) = (3.13) 

n=0 

To ensure that Se C S^^e, we just need to choose a parametrisation of S such that 0, 0) = 
e(t) which fixes the Taylor coefficients S^'^''^' ''^\0,...,0) = e^"^\0). For the finite number 
k = \E{v) \ — 1 of remaining edges at v, we can now use the freedom in choosing the other Taylor 
coefficients to assure that there are no (beginning segments of) other edges contained in Sy^e 
[24:\. In particular, only a finite number of Taylor coefficients is involved. 

Now we state that the intersection properties of a finite number of transversal edges at v 
with any (sufficiently small) surface S are already fixed by a finite number of Taylor coefficients 
of S. We will discuss the case D = 3 for simplicity, higher dimensions are treated analogously. 
Locally around v we may always choose coordinates such that the surface is given by z = 0, 
S{x,y) = (x,y,0). The edge e contained in the surface is given by e{t) = {x{t),y(t),0) and for 
any transversal edge at v we find e'{t) = {x' (t) , y' (t) , z'{t)) where z'{t) = ^^"_^-^, z'("~-^)(O) + 0(P), 
and n < oo since otherwise e' would be contained in S. The sign of the lowest non-vanishing 
Taylor coefficient z '^"■"^^(0) determines if the edge is "up"- or "down"-type locally. Set = 
T^8L^e'eE{v),e'^e ("■)> ^^icl obviously N < OO. Thus, we can e.g. by modifying 5(^''^)(0, 0) choose 
the surface S{x,y) = (x,y,ibx^), which locally has the same intersection properties with the 
edges e' G E{v), e' ^ e and certainly does not contain e anymore. 

Coming back to the general case considered before, there always exists N < oo such that we 
can change £'(^''^'■••''^''(0, ...,0) without modifying the intersection properties of any of the edges 
e' G E{v),e' / e, in particular the "up"- or "down"-type properties are unaffected. However, 
the edge e no longer is of the inside type, but becomes either "up" or "down" (depending on 
whether 5(^'°'-'°)(0, 0) is scaled up or down and on the orientation of S). In general, new 
intersection points v' G E{v) Ci S,v' ^ v may occur when modifying the surface in the above 
described way, but we may always make S smaller to avoid them. 

Now choose a pair of surfaces 5, S for the edge e such that it is once "up"- and once 
"down" -type to obtain the desired result 

Yij{S) - Yjj{S)] p;f^ = 2p;R<jjf^. (3.14) 

Case (h): In the case that there is a partner e which is a analytic continuation of e through 
V, we cannot construct an analytic surface (without boundary) Sy^e containing a beginning 
segment of e and not containing a segment of e. However, we can construct an analytic surface 
Sv,{e,e} containing (beginning segments of) e, e and sharing the remaining properties with S^^e 
above. The method is the same as in case (a) |24| . Again, there always exists N < oo such 
that we can change S^^''^''"'^\0, 0) without modifying the intersection properties of any of the 
edges e' G E{v),e' / {e,e}, and such that both edges e, e become either "up" or "down"-type. 
Moreover, if we choose A^ even, then e, e will be of the same type with respect to the modified 
surface, while for A^ odd one edge will be "up" and its partner will be "down". Calling the 



modified surface S for A^ even and S for A^ odd, we find with the same calculation (3.14) as in 
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case (a) the desired result. 

This furnishes the proof of the above lemmsj^ □ 
Choosing the surfaces as described above, we find that the fohowing hnear combination 

^ (C^{S, S', x) - C^{S, S', x) - C^{S, S', x) + C^{S, S', x)) p*^f^ 



.IJKLMj^e^j^e^j^ (3.15) 



4 



proves the necessity of the building blocks. Using the fact that the edge representations are 
already simple, we can rewrite the building blocks as 



Tie rye' f _ ^ \ ( jye , T?e' \( ne . T>e' \ pe pe pe' jje' 



f' 
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2 

= 2^^iJ + ^[ij)(^KL] + ^KL])fi -■ 2^fjKd-i- (3-16) 

We proceed by showing that the building blocks are anomalous, starting with the case D = 3. 
We calculate for e 7^ e' / e" 7^ e 

^fjKL^^'^^'^^ ^ABCD ^^fFK{Re")AB{ReY'^{Re')^C^ (3.17) 

where we used the notation := n\ 5^j^6j\..Sj"y To show that this expression can not 



be rewritten as a linear combination of the of building blocks (3.16), we antisymmetrise the 
indices [ABU], [ABKC] and [IJKC] and find in each case that the result is zero. Therefore, 
a simplicity building block can not be contained in any linear combination of terms of the type 



(3.17). For D > 3, we have 



JJKLM /.ee' ABODE /.e'e" ^ABCE(n \ (t>\IJ(v \K 1 8^ 

^IJKL^^ ^ABCD\ ^ (>jjKM^^<^")AB[ne) [He') C- (3.18) 

Choosing M = E fixed, the anomaly is the same as above. A short remark concerning the 
terminology "anomaly" is in order at this place. Normally, the term anomaly denotes that a 
certain classical structure, e.g. the constraint algebra, is not preserved at the quantum level, e.g. 
by factor ordering ambiguities. The non-commutativity of the simplicity constraints however can 
already be seen at the classical level when using holonomies and fiuxes as basic variables. Thus, 
one could argue that it would be more precise to talk of a quantisation of second class constraints. 
On the other hand, since the holonomy-flux algebra is an integral part of the quantum theory 
and at the classical level it would be perfectly fine to use a non-singular smearing, we will 
nevertheless use the term anomaly to describe this phenomenon. 

Independently of the terminology chosen, we cannot quantise the simplicity constraints acting 



on vertices using the Dirac procedure since this will lead to the additional constraints (3.18) 
being imposed. The unique solution to these constraints has been worked out in |23j and is 
given by the Barrett-Crane intertwiner in four dimensions and a higher-dimensional analogue 
thereof. Several options are at our disposal at this point. Looking back at our companion 
paper [2], one could try to gauge unfix this second class system to obtain a first class system 
subject to only a subset of the vertex simplicity constraints. In this process, one would have 
to pick out a first class subset of the simplicity constraints which has a closing algebra with 
the remaining constraints. The construction of a possible choice of such a subset is discussed 



^This also establishes that the right invariant vector fields i?f j are not only contained in the Lie algebra 
generated by the flux vector fields Y{S), but are already contained in the flux vector space, which to the best of 
our knowledge has not been shown. 
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in our companion paper [25]. While the proposed subset is first class with respect to the other 
constraints, it suffers from the fact that the choice is based on a certain recoupling scheme and 
that a different choice of the recoupling scheme results in a different first class subset. This 
is not a problem for the theory itself, but it seems problematic when constructing a unitary 
map to SU(2) spin networks in four dimensions, as discussed in [25]. Another possibility is the 
construction of a Dirac bracket, which however would result in a non-commuting connection 
and the non-applicability of the LQG quantisation methods. The use of a weak implementation 
in the sense of Gupta and Bleuler is discussed in our companion paper [25j. While the results 
obtained in the context of the EPRL spin foam model can be also used in the canonical theory 
(up to certain subtleties discussed in [25]), they rely on specific properties of SO (4) which do 
not extend to higher dimensions. 

While equivalent at the classical level, the master constraint introduced in [26] allows to 
quantise also second class constraints by a strong operator equation. Due to the second class 
nature, one expects the master constraint operator to have an empty kernel or at least a kernel 
which is too small to describe the physical Hilbert space. Since we know that the Barrett- 
Crane intertwiner is a solution to the strong imposition of all vertex simplicity constraints, we 
are in the second case. In order to find a larger kernel of the master constraint, one modifies 
it by adding terms to it which vanish in the classical limit, i.e. performs ^-corrections. The 
merits of this procedure are exemplified by the construction of the EPRL intertwiner fSj in four 
dimensions, which results from a master constraint for the linear simplicity constraint upon 
/i-corrections. A simplification arising in the treatment of the linear simplicity constraints, see 
e.g. [8] our companion paper [25], is that they act individually on every edge connected to the 
intertwiner. On the other hand, the quadratic constraints act on pairs of edges and the resulting 
algebraic structure of the master constraints is thus very different. Since we are not aware of a 
suitable solution for the quadratic vertex master simplicity constraint, we will contend ourselves 
by giving a definition of this constraint operator. The task remaining for solving the vertex 
simplicity master constraint operator is thus to find a proper ^-correction which results in a 
physical Hilbert space with the desired properties, e.g. that there exists a unitary map to SU(2) 
spin networks in four dimensions. 

A general simplicity Master constraint is given by 

Mn*f -„* V „e"e"'MNOPA,ee' Ae"e"' r i q^ 

— C(,ei IJKL ^IJKL^MNOPJl [6.L^) 

e,e',e",e"'&E{v) 

with a positive matrix Cg^f Yjkl^ ^ which we will choose diagonal for simplicity, c^gf Yjkl^ ~ 
jfCee'i^^^g <5g/) <5/j;^i'^- The diagonal elements c^e' can be chosen symmetric because of the sym- 
metry of the building blocks. We choose c^e' = lVe,e',e / e' and Cge = since the edge 
representations are already simple, leading to the final version of the Master constraint we 
propose, 

M,p;f, = p; Yl ^fjKL^fjKLf,- (3.20) 

Cylindrical consistency and essential self-adjointness follows analogously to the case of C{S, S' , x) 
in section fS. 2. 1[ 

For the case of S0(4), we can use the decomposition in self-dual and anti-selfdual generators 
to find that e^-^^^R'ijRi^L = ■ - Jl ■ J!_', which implies 

IJKL A ee' / re i je' 1 I \ re 



e'^^^Afj^i = + • + jn - + • + Jl =: - A!! . (3.21) 
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This leads to the Master constraint 

M,p;f^=p; (Af -2Af A"=' + A!_^'A1^')/^, (3.22) 

where + and — now label independent copies of S0(3). Thus, we can calculate the matrix ele- 
ments of this constraint in a recoupling basis analogously to the standard LQG volume operator 
matrix elements [27j . 

As mentioned before, alternative routes to deal with the vertex simplicity constraints will be 
the subject of [25]. 

3.3 Diffeomorphism Constraint 

The diffeomorphism constraint can again be treated in exact agreement with the (3 + 1)- 
dimensional case. To solve the diffeomorphism constraint, one proceeds as follows. Consider the 
set of smooth cylindrical functions $ := CyY^ (A/G) which can be shown to be dense in T-L^. By 
a distribution -(/^ S on $ we simply mean a linear functional on $. The group average of a 
spin-network state ^ ^ is defined by the following well-defined distribution on <1> 

^h]A,c--= E <^7'Ac'->' (3-23) 

7' 6 [7] 

where [7] denotes the orbit of 7 under smooth diffeomorphisms of a which preserve the analyticity 
of 7 including an average over the graph symmetry group (see, e.g., [28j for technical details). 
Since we already solved the simplicity constraint on single edges, we can restrict attention to 
spin network states with edges labelled by simple SO{D + 1) representations, Ae = {Ne,0, ■■■)■ 
The group average [/] of a general cylindrical function / is defined by demanding linearity of 
the averaging procedure, i.e. first decompose / into spin-network states and then average each 
of the spin-network states separately. An inner product for the diffeomorphism invariant Hilbert 
space can be constructed. We will not give details and refer the reader to |16l I28j . 

4 Geometrical Operators 

4.1 The D — 1 Area Operator 

The area operator was first considered in and defined mathematically rigorously in the LQG 
representation in [30]. In [1], the results of [30] are generalised for arbitrary dimension D. Using 
the classical identity we can basically copy the treatment found there. Let 

5 be a surface and X : Uq ^ S the associated embedding, where Uq is an open submanifold of 
]^D-i^ Then the area functional is given by 

Ar[S] := / d^-^uyJdet{[X*q] {u)). (4.1) 

Introduce Uq = Ujj^uU, a partition of Uq by closed sets U with open interior, U being the 
collection of these sets. Then the area functional can be written as the limit as |?7| — )• 00 of the 
Riemann sum 

At[S] := Jl7rjj{Su)7r^-'{Su), (4.2) 
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where Sjj = X{U) and t^ij{Su) is the electric flux with choice n^'^ = S^j^dj^-^, which has been 

quantised already. Let / G Cyl^(^), choose a representative and, using the known action of 
the quantised electric fluxes, obtain as in the (3 + 1) -dimensional case 



x£{enS;e£E{^s)} 



\ 



e£-B(7s),a;6<9e 



7' 



(4.3) 



where 75 ;^ 7 is an adapted graph. The family of operators Ar^[5] has dense domain Cyl^(^). 
Its independence of the adapted graph follows from that of the electric fluxes. Moreover, the 
properties of the area operator like cylindrical consistency, essential self-adjointness and dis- 
creteness of the spectrum can be shown analogously to [1]. 

The complete spectrum can be derived using the standard methods. We use 



^ ee-E(7s),a;e9e 



D ^ j-tX,up\2 „ / jjX,down\'^ ( j-,x,up j-)X,down\' 

^ K^ij ) + ^ y^u ) - [^ij + ^ij ) 



A itp \down _j_ \ \ up+down 

2 



(4.4) 



where the As are mutually commuting primitive Casimir operators of SO{D + 1). Thus their 
spectrum is given by the Eigenvalues Ajr > 0. We have to distinguish the cases D + 1 = 2n even, 
N 9 n > 2 and D + 1 = 2n -|- 1 odd, n G N. In a representation of SO(D + 1) with highest weight 
A = (ni, ...,nn), fii G No, we find for the eigenvalues of the CasimiJ^ A := — ^A/jX^"^ 



AfA := Atta'Ua 



AfA := K^v^ 



j=2 i<j 



i=l 



j=2 i<j 



v/^ for S0(2n), 

VA forSO(2n + l), 



i=l 



1=1 



(4.5) 



where we used the following notation 

n-2 



E"j + 



]=i 



n-n-1 + n-r, 
2 



i<{n- 2); fn-i = ; fn = for S0(2n), 



fi = "^rij + i < {n - I)] fn 



for S0(2n + 1), 



(4.6) 



such that /i > /2 > ... > fn- Note that the above formulas hold for general irreducible 
Spin(Z)+l) representations. Irreducible representations of S0(Z)+1) are found by the restriction 
that all fi be integers. Denoting by 11 a collection of representatives of irreducible representations 
of SO(D -|- 1), one for each equivalence class, we find for the area spectrum 

Spec(Ar[S]) = | ^ E ^2A,i + 2A^^ " ; G N, vr^ , vr^ , ^ n, G 4 ® • (4-7) 

*Note that R''' = 1/2X^^ , such that X^'^ fulfil the standard Lie algebra relations without the factor 1/2 
appearing in ( 2.8 1. 
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Note that the above formulas (4.5) significantly simplify if we restrict to simple representations, 
Ao = (A^,0,0,...), 



AvAo = iV(iV + 2n - 2)uAo =N{N + D -l)vA, for S0(2n), 

AvAo = N{N + 2n + l-2)vAo=N{N + D -l)vA^ forSO(2n + l). (4.8) 

We cannot use this simplified expression for the SO{D + 1) Casimir operator in the general case 



(4.7), since in the decomposition of a tensor product of irreducible simple representations usually 
non-simple representations will appeai|^ but we can use it for a single edge. When acting on a 
single edge e = ei o (62)^^ intersecting S transversally, we know that due to gauge invariance 

{Rfj - RfjY K = A {RfjY he = -2N{N + D- l)he. (4.9) 



The action of the area operator on a single edge e, e n 7^ is thus given by 

D-l 

) 



APe[S]plhe = Kh{3^/N{N + D-l)plhe = 167r/3 ^ X ^/N{N + D-l)plhe, (4.10) 



where lp^~^^^ := ^ y ^'^^J;^^' is the unique length in D + 1 dimensions, and k = 16ttG^^^^^ /c^ 

in any dimension, where G^^^^^ denotes the gravitational constant. Note that for D = 3, we 
find the factor y^N{N + 2) in the area spectrum of an edge stemming from irreducible simple 
representations of SO (4). Replace the non- negative integer labelling the weight by = 2j, j 
half integer, to find the factor 2y^ j{j + 1) of SO (4) spin foam models, which coincides with the 
usual spacing in (3 + l)-dimensional LQG, 

APe[S]plhe = 2Khpy%Tl)p:K = 327r/3 (liF+^'^Y'^ X Vi(i + 1)p:K. (4.11) 

In standard LQG, instead of the gauge group SO (3) one extends to the double cover Spin(3) = 
SU(2) and allows also for half integer representations. Note that in our case, we cannot allow 
for general Spin(Z) + 1) representations at the edges, since the edge simplicity constraint is not 
satisfied in representations of Spin(L' + 1) which are not as well representations of SO{D + 1), 
L> > 3 1231. 



4.2 The Volume Operator 

The derivation of the volume operator is analogous to the treatment in [1] and requires only a 
slight adjustment. 

The volume of a region R is classically measured by 

ViR) := [ d^x^g, (4.12) 
JR 

where y/q has to be expressed in terms of the canonical variables. The derivation is performed 
for /3 = 1, the general result is obtained by multiplying the resulting operator by f3^^^^~^\ 

^For the tensor product of two irreducible simple representations of SO(n) holds [311 132j (w.l.o.g. M > N) 
[M, 0, .., 0] ® [TV, 0, .., 0] = Ek=o El=o'^[M + N-2K-L,L,0,..,0]. 
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4.2.1 D + 1 Even 



Let n = (D—l) /2. Let xa{Pi x) be the characteristic function in the coordinate x of a hypercube 
with centre p spanned by the D vectors A* := A*n*, i = 1, . . . , D, where n* is a normal vector in 
the frame under consideration and which has coordinate volume vol = . . . A^ det(n^, . . . , n^) 
(we assume the vectors to be right-oriented). In other words, 

XA{p,x) = l[Q(^-\<n\x-p>\] (4.13) 

i=i ^ ^ 

where < •, • > is the standard Euclidean inner product and Q{y) = 1 for y > and zero otherwise. 
We will use lower indices (A], . . . , Af ) to label different hypercubes. It will turn out to be 
convenient to label the D edges appearing in the following formulae by e, ei, . . . , e„, e'l, . . . , e^. 
We consider the smeared quantity 

7r(]3, Ai, . . .,Ad) 



(7 



V0l(Ai)...V0l(Ai5) J a 

XAi (p, xi)xA2 (2p, X1+X2)... XAo {Dp, xi + ... + xd) 

^eaai6i...a„6„e7J/i Ji/2J2.../nJ„7r"^''7r"^'i^%''i-'iKi . . . tt""^"^" A-^"x„. (4.14) 
Then it is easy to see that the classical identity 

V{R)= lim ... lim f d^p Mp, Ai, . . . , Ad)\'^ (4.15) 

holds. The canonical brackets 

{Aauix), 7r^^^(2/)} = 25^ (x - y)Slsf6^^ (4.16) 
give rise to the operator representation 

I SAbKL 

while the connection acts by multiplication. 

Let a graph 7 be given. In order to simplify the notation, we subdivide each edge e with 
endpoints v,v' which are vertices of 7 into two segments s, s' where e = s o (s')^^ and s has an 
orientation such that it is outgoing at v'. This introduces new vertices sH s' which we will call 
pseudo-vertices because they are not points of non-semianalyticity of the graph. Let £'(7) be 
the set of these segments of 7 but ^^(7) the set of true (as opposed to pseudo) vertices of 7. Let 
us now evaluate the action of 

7r»"(p, A) := ^ d^x x{p, ^)^"" (4-18) 

on a function / = cylindrical with respect to 7. We find (e : [0, 1] — >• a, t — )■ e{t) being a 
parametrisation of the edge e) 

7r«^^(p, A)/ = Yl ^^A(^>,e(^))e«(^)tr (^[/,,(0,t)r^-^^,(i, 1)]^ (4.19) 

e£B(7) 
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Here we have used (1) the fact that a cyhndrical function is aheady determined by its values on 
A/Q rather than A/Q so that it makes sense to take the functional derivative, (2) the definition 
of the holonomy as the path-ordered exponential of A with the smallest parameter value to 
the left, (3) A = dx'^AaijT^'^ where r^"^ G so{D + 1) and we have defined (4) ti{h^d/dg) = 
hAsd / Bab, A, B,C, . . . being SO{D + 1) indices. The state that appears on the right-hand side 



of (4.19) is actually well-defined, in the sense of functions of connections, only when A is smooth 



for otherwise the integral over t does not exist, see [33] for details. However, as announced, we 
will be interested only in quantities constructed from operators of the form 
the limit of shrinking A — )• to a point has a meaning in the sense of Ti 



^2{A/Q,dfio) and 

therefore will not be concerned with the actual range of the operator (4.19) for the moment. 

We now wish to evaluate the whole operator Tr{p, A^, . . . , A-^) on /. It is clear that we obtain 
D types of terms, the first type comes from all three functional derivatives acting on / only, the 
second type comes from D — 1 functional derivatives acting on / and the remaining one acting 



on the trace appearing in (4.19), and so forth 



The first term (type) is explicitly given by 
^(p,Ai,...,Ad)/ 

1 (ih)^ 



(4.20) 



2L>!vol(Ai)...vol(Ai5 

XAi(p,2;i)xA2(2p,2;i + X2 

tl([he{0,ty'^he{tA)f 



^aaibi.. (^IJhJll2J2-InJn 



■ XAd (Dp, xi + 



'[0,1P 

+ XD)e''{t)e1'{h] 



dt dti . . . dtn dt'i . . . dt[ 



E 

ei,...,eD&E(^) 
^n"(^n)e'i^^(tl) . . . e'J'"{tn) 



d 



dheiOA) 



T 



tr (^[/lei(0,tl)r^l^l/lei(tl,l)] 
tr f[/le„(0,t„)r^"^"/le„(in,l)]'^ 



d 



dheAOA] 
d 

dhZiOA] 



tr 



tr 



d 



a^(o,i) 

T d 



dKAQA] 



The other terms are vanishing due to either the same symmetry / anti-symmetry properties as 
in the usual treatment or the simplicity constraint in case the first derivative is involved. 
Given a D-tuple ei . . . of (not necessarily distinct) edges of 7, consider the functions 



Xex,...,eD{h-, ■■■■.to) ■■= ei(ti) + . . . + eoitD)- 
This function has the interesting property that the Jacobian is given by 



det 



X 



D 



)iti,...,tD) 



,tD) 



':ai...aD 



hit 



\0-D 



(4.21) 



(4.22) 



which is precisely the form of the factor which enters the integral (5.8) 
We now consider the limit A^, . . 



0. The idea is that all quantities in (5.8) are 



meaningful in the sense of functions on smooth connections and thus limits of functions as 
A — )• are to be understood with respect to any Sobolev topology. The miracle is that the final 
function is again cylindrical and thus the operator that results in the limit has an extension to 

all of A/g. 
Lemma. 

For each D-tuple of edges ei, . . . , e^i there exists a choice of vectors n\, . . . , n\),n1, . 
a way to guide the limit A\, A2, . . . , A^ — )• such that 



, and 



det 



d{ti, ...,tD) 



XAi ip,ei)... XAd {Dp, ei + ... eD)6e^,. 



(4.23) 
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vanishes if 



(a) if ei, . . . ,eo do not all intersect p or 

= (which is a diffeomorphism invariant statement). 



(b) det 1^ 



d{ti,...,tD) 



Otherwise it tends to 



D 



(4.24) 



i=l 



Here we have denoted by Oei,...,eD (p) the trace(s) involved in the various terms of (5.8) 



We conclude that (5.8) reduces to 



lim 7r(p,Ai,...,AD)/ 



E 



ei,...,ei3 



{ih)^s{ei, . . . ,eg) 
2^D!vol(Ai)...vol(Ai5-i 



-XAi (P, v)... XAn-i (P> l')Oei,...,eo (0, . . . , 0), (4.25) 



where v on the right-hand side is the intersection point of the -D-tuple of edges and it is under- 
stood that we only sum over such D-tuples of edges which are incident at a common vertex and 
s{ei, . . . , eo) '■= sgn(det(ei(0), . . . , e£)(0))). Moreover, 



O. 



ei,...,e£) 



(0, . . . , 0) = ^eiji,j,i,j,..j„j„Ri'Rif^R-'^}K^ ■ ■ . Rf^-KlK. 



and 



R. 



ij ._ 



R'-\he{0,l)) ■.= tl({T''he{0,l)) 



d 



5/ie(o,i: 



(4.26) 



(4.27) 



is a right-invariant vector field in the r^"^ direction of SO{D + 1), that is, R{hg) = R{h). We 
have also extended the values of the sign function to include 0, which takes care of the possibility 
that one has D-tuples of edges with linearly dependent tangents. 

The final step is choosing Ai = . . . = /S.d-i and exponentiating the modulus by 1/{D — 1). 
We replace the sum over all D-tuples incident at a common vertex by a sum over all 

vertices followed by a sum over all -D-tuples incident at the same vertex Ylv<^v('y) Sein r\eo=v 
Now, for small enough A and given p, at most one vertex contributes, that is, at most one 
of Xa{v-,p) 7^ because all vertices have finite separation. Then we can take the relevant 
Xa(P; v) = xa{p, v)'^ out of the exponential and take the limit, which results in 



V{R) 
V{p) 

Qei,...,er> 



[ d^p|det'(^)|7= / d'^pVip)^, 
Jr Jr 



D 



i^ 

D\ 



X] s(ei,...,ez))gei,...,eo 

ei,...,e£)£-E(7), ein...neo=D 



1 



'-,eijHJM...i^j„Ri'Rif'R'jK, ■ . . Kl^'-RilK^. 



1 

D-1 



(4.28) 
(4.29) 

(4.30) 
(4.31) 
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4.2.2 D + 1 Odd 



The case D + 1 uneven works analogously, except that the expression for det(q) is changed a 
bit. With n = D/2, the result is 



V{R) 
V{p) 



D,7 



?ei,...,eo 



rf>|det(g)(p)|^= / d^'pVip)^, 
R JR 



D 



vev{-f) 



D\ 



ei,...,e£,&E{-y), e\n...neo=v 

1 

yi y \2D-2 



eiiMJ,...i^j„Rif'Ri}K, . . . Rif-R'/K^. 



(4.32) 
(4.33) 

(4.34) 

(4.35) 
(4.36) 



4.2.3 More Results and Open Questions 

The derivations of cylindrical consistency, symmetry, positivity, self-adjointness and anomaly- 
freeness given in [3] generalise immediately to the higher dimensional volume operator. The 
question of uniqueness of the prefactor [311 ESj in front of the expression under the square root 
of the volume operator or the computation of the matrix elements [36l [Ml EZl EHl |39] have not 
been addressed so far, however these are not necessary steps in order to use the volume operator 
for a consistent quantisation of the Hamiltonian constraint in what follows. We leave these open 
questions for future research. 



5 Implementation of the Hamiltonian Constraint 
5.1 Introductory Remarks 

The implementation of the Hamiltonian constraint will follow along the lines of see [40] for 
original literature and details. In our companion papers [Tl|2], we derived the classical expression 

n:=^^ (-^P-^Ue + fj, (/^)Z)^^ - + 1) K^,K,jE^y^E'\'^'^ , (5.1) 

where a,b,c, . . . are spatial indices and I, J, K, . . . are so{D + 1) indices. In order to have a well 
defined quantum version of this constraint, we have to express it in terms of holonomy and flux 
variables. As in the 3 + 1-dimensional case, the volume operator turns out to be a cornerstone 
of the quantisation. 

At first, we will introduce a graph adapted triangulation of a in order to regularise the Hamil- 
tonian constraint. Next, classical identities to express the Hamiltonian constraint in terms of 
holonomies and fluxes are derived. Since the complete expression for the Hamiltonian constraint 
will turn out to be rather laborious to write down, we will derive the regularisation piece by 
piece. Next, we show how to assemble the regularised pieces to the complete constraint and 
describe the quantisation. Finally, we construct a Hamiltonian Master constraint in order to 
avoid some of the usual difficulties associated with quantisation. 
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5.2 Triangulation 



A natural choice for a triangulation turns out to be the following (we simplify the presentation 
drastically, the details can be found in |40j): given a graph 7 one constructs a triangulation 
T(7, e) of a adapted to 7 which satisfies the following basic requirements. 

(a) The graph 7 is embedded in T(7, e) for all e > 0. 

(b) The valence of each vertex v of 7, viewed as a vertex of the infinite graph r(7, e), remains 
constant and is equal to the valence of v, viewed as a vertex of 7, for each e > 0. 

(c) Choose a system of semianalyticj^ arcs ^ g g/, one for each pair of edges e,e' of 7 incident 
at a vertex of 7, which do not intersect 7 except in its endpoints where they intersect 
transversally. These endpoints are interior points of e, e' and are those vertices of T(-y, e) 
contained in e, e' closest to v for each e > (i.e., no others are in between). For each e, e' > 
the arcs a^^ge'' ^^vee' diffeomorphic with respect to semianalytic diffeomorphisms. 
The segments e, e' incident at v with outgoing orientation that are determined by the 
endpoints of the arc ^ g g/ will be denoted by ,^ g, ^ g, respectively. Finally, if (/> is a 
semianalytic diffeomorphism then s^^^^^^^^^^^^^y «^(7),</,(,.),<^(e),0(e') and Hsy,v,e), Ha'^,v,e,e') 
are semianalytically diffeomorphic. 

(d) Choose a system of mutually disjoint neighbourhoods U^^^, one for each vertex v of 7, and 
require that for each e > the at,,^„, are contained in U^„. These neighbourhoods are 

nested in the sense that [/^ j, C U^^^ if e < e'. and limg_!.o U^^^ = {v}. 

(e) Triangulate U^^^ by D-simplices A(7, ei, . . . , e/)), one for each ordered i^-tuple of distinct 
edges ei, . . . incident at bounded by the segments ,^ g^, . . . , s^^^^g^ and the arcs 
«7,,;,ei,e2'«7,t;,ei,e3>--->«7,«,eo-i,eo {D{D-l)/2 arcs) from which loops a^;„;ei,g2, etc. are 
built and triangulate the rest of a arbitrarily. The ordered D-tuple ei, . . . , e^) is such that 
their tangents at v, in this sequence, form a matrix of positive determinant. 

Requirement (a) prevents the action of the Hamiltonian constraint operator from being 
trivial. Requirement (b) guarantees that the regulated operator H'^{N) is densely defined for 
each e. Requirements (c), (d) and (e) specify the triangulation in the neighbourhood of each 
vertex of 7 and leave it unspecified outside of them. 

The reason why those D-simplices lying outside the neighbourhoods of the vertices described 
above are irrelevant will rest crucially on the choice of ordering with [hj^, V] on the rightmost: 
if / is a cylindrical function over 7 and s has support outside the neighbourhood of any vertex of 
7, then V{'jU s) — V{'y) consists of planar at most four-valent vertices only so that [hj^, V]f = 0. 

We will define our operator on functions cylindrical over coloured graphs, that is, we define 
it on spin network functions. The domain for the operator that we will choose is a finite 
linear combination of spin-network functions, hence this defines the operator uniquely as a 
linear operator. Any operator automatically becomes consistent if one defines it on a basis, the 
consistency condition simply drops out. 

The volume operator will appear in every term of the regulated Hamiltonian constraint. 
We will choose a factor ordering such that the Hamiltonian constraint acts only on vertices. 
It is therefore sufficient to regularise the constraint at vertices. As in the usual treatment, we 
use the tangents to the edges at a vertex as tangent vectors spanning the tangent space of 
the spatial coordinates. To emphasise this, we will abuse the notation in the following way: 
Let ea(A) denote the D edges incident at the vertex v of an analytic L'-simplex A G T(7,e). 

^Semianalyticity is a more precise version of piecewise analytic. See [17] for complete definitions. 
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The matrix consisting of the tangents of the edges ei(A), . . . , eD(A) at v (in that sequence) 
has non-negative determinant, which induces an orientation of A. Furthermore, let aab be 
the arc on the boundary of A connecting the endpoints of ea(A), ef,(A) such that the loop 
aab{^) = ea(A) o aabi^) o e;,(A)^^ has positive orientation in the induced orientation of the 
boundary for a < 6 (modulo cyclic permutation) and negative in the remaining cases. 

5.3 Key Cleissical Identities 

The following classical identities are key for the rest of the discussion. 

5.3.1 i:> + 1 > 3 Arbitrary 

We observe that 

VQ^aiAx) := -{D - l){Aaij, V{x, e)}, (5.2) 

where V{x,e) := / d^y x<i{x,y)y/q is the volume of the region defined by Xe{x,y) = 1 measured 
by Qab and Xe{x, y) = 11^=1 0(e/2 — — ?/"|) is the characteristic function of a cube of coordinate 
volume with centre x. Also, 

n\x)nj{x) « (7r«^'(x)7r„;^j(x) - j) . (5.3) 

We can write the KKEE terms in the same way as in the usual 3 + 1-dimensional case, 
using 



Further, 



K{x) := E''\x)Kai{x) ?a {UEix), V{x, e)} . (5.4) 



E'^\x)Kai{x) « ^^^T,b^L{x) {AaKL{x), {Hfi[l](x, e), V{x, e)}} (5.5) 



gives us access to all the needed terms. 

5.3.2 D + 1 Even 

Let n= {D — l)/2. It is easy to see that 

''""(^) ^^^^e"''i^i-^"'^"e-^-^-^i-^i-^"^"sgn(dete)(x) 

nihK^ix)Trc^Ji^^{x) . . .7rb„/„K„(ar)7rc„j„^"(x)^^"^(x). (5.6) 

The sign of the determinant of where the internal space is the subspace perpendicular to 
is accessible through 

Sgn(det(ei))(x) « Ji-Jn^n,aai6i...a„6„^D-l^^^^(^) 

vrai/iKi(a;)7rfeij,^i(a;) • ■ ■ T^a„inK„ix)7rb^j^^" (x). (5.7) 
For the Euclidean part of the Hamiltonian constraint, we need 

~77?i -e"^^"i^i-""-i^"-ie-^-^™-^i--^"-i-^«-isgn(det e)(x) (5.8) 

y/Q 4:[D — 2). 

T^cKLix)7raihKi ix)7rb^j^^^ (x) . ..Tran-iI„-iKn-iix)7rbn-iJn-l^"~Kx)VQ^~^{x). 
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Regarding quantisation, wc have to choose a classical expression for ^ X i^)- The above 
expression would be favourable by arguments of simplicity if it would not contain the additional 
factor of sgn(det(e^))(a;) which has to be accounted for. Therefore, we can equally well express 
the two factors of tt"^"^ separately and absorb the inverse square root into volume operators. 

5.3.3 D + 1 Odd 

Let n = {D — 2)/2. With only minor modifications of the D + 1 even case, we get 



T^bihK,{x)7rc,j,^'{x) . . .nb„i„Kn{x)7rcnJn^"{x)Vq^ \x) (5.9) 



with 



n^{x) « -^e"i''i-""+i^"+ie^^i'^i--^"+i-^"+isgn(dete)(x)^^"^(x) 

T^aihKi{x)Trb,Ji'^'{x) . . .7r„„+i/„+iif„+i(x)7r6„+,j„+l-^"+l(x). (5.10) 

For the Euclidean part of the Hamiltonian constraint, we need 

^ ~ 2(L>-2)!' ' sgn^detej 

riKT^aihR^T^biJi^^ ■ ■ ■ 7ra„7„Kn7r6„j„^"^^~^ (5.11) 

and observe that the factor of sgn(det(e^))(x) is canceled by another such factor coming from . 
The Euclidean part of the Hamiltonian constraint therefore has the same amount of complexity, 
measured by the "number of involved operators" , in even and odd dimensions. 

5.4 General Scheme 

The basic idea of the regularisation of the Hamiltonian constraint operator is to approximate 
the constraint operator on the graph adapted triangulation and then to take the limit of an 
infinitely refined triangulation. For this procedure to work, it is mandatory that the constraint 
operator has a density weight of +1. A typical term of the classical Hamiltonian constraint (or 
any other operator one wants to regulate) will, after using the above classical identities, consist 
of 

• an integral d^x, 

• n G No spatial e symbols, 

• factors of Aaij{x), 

• Poisson brackets involving a factor of Aaij{x) as one of its two arguments as well as either 
the volume of a neighbourhood of x, the Euclidean part of the Hamiltonian constraint 
smeared with unit lapse over a region containing x, or the Poisson bracket of the Euclidean 
part of the Hamiltonian constraint with the volume, smeared as before, as the other 

argument, 

• field strength tensors, 

• a factor of ■s/q^~"', 
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• (covariant) derivatives. 

Operators that are well defined on the kinematical Hilbert space are holonomies and the volume 
operator. We will show in the following that we can construct the Euclidean part of the Hamilto- 
nian constraint operator, which gives us access to the remaining part of the constraint operator. 
As a start, it is therefore mandatory to write the Euclidean part of the Hamiltonian constraint 
in terms of holonomies and volume operators. We stress that we do not quantise the tt""^^ as 
flux operators, which would also be possible. The reason is that the Hamiltonian constraint 
operator would not simplify significantly by using fiuxes instead of derived flux operators. On 
the other hand, the appearance of fluxes only through volume operators can be seen as a certain 
simplification. Anyhow, different regularisations are possible and the discrimination between 
different regularisations has to be considered in the semiclassical limit. 

We begin with rewriting the integral. Given a -D-tuple of edges (ei, . . . ,6/)) incident at v 
with outgoing orientation consider the D-simplex A'^(7, ei, . . . , e^) bounded by the D segments 
^,i;,ei' • • ■ ) ^'^^,v,eo incident at v and the D{D — l)/2 arcs al^^^^ea^e,,-! ^ < a < h < D. We now 



define the "mirror images" 



'^■y,v,p,p' (^) 



a. 



"^v - 4,„,p(t), 

y,v,p,p' (0 ' 



2v 



•,v,p,p' (0 



2t\v 



/(I)], 



a7,„,p,p'(i) + 2tb-4,„y(l)], 



(5.12) 



where p ^ p' G ei,...,eD and we have chosen some parametrisation of segments and arcs. 
Using the data (5.12) we build 2^ — 1 more "virtual" D-simplices bounded by these quantities 
so that we obtain altogether 2^ D-simplices that saturate v and triangulate a neighbourhood 
Let J, be the union of these neighbourhoods as we vary the ordered -D-tuple 



of edges of 7 incident at v. The U^^^, v G ^(7) were chosen to be mutually disjoint in point (d) 
above. Let now 



U. 



7,i;,ei,...,e£, 



"-^7,?;,ei,...,e/3 5 

^- U 

?^ey(7) 



(5.13) 



then we may write any classical integral (symbolically) as 



^7 t,ev(7) ^ ^ 



c/^^ ^ E{v) 
„ev(7) ^ ^ 



ti=fe(ei)n...nb(e£3) 
-i;=fe(ei)n...n6{e£,) 



+ 



^7,u,ej,...,e25 



'-i,v,ei,...,eQ , 



7,f ■ei,...,e^ 



7,i',ei,...,ej5 



(5.14) 



where in the last step we have noticed that classically the integral over K^^^^ei,...,ei3 converges 
to 2^ times the integral over ^^^v,ei,...,eo' ~ means approximately and E{v) = ("jj"^) with 
n{v) being the valence of the vertex. Now when triangulating the regions of the integrals over 



^,v,ei,...,eo 



and in (5.14), regularisation and quantisation gives operators that vanish on 



because the corresponding regions do not contain a non-planar vertex of 7. 
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As a next step, we approximate the integral 



J^^ d^x9{x)f^^e''9iv) (5.15) 



"f,v,ei , 



for some function g{x). Here we assumed the coordinate length of each segment s^, ,,, to be 
e. The general case of arbitrary coordinate length works analogously, since the factors of e will 
be hidden in holonomies and derivatives contracted with an epsilon symbol which addresses 
each segment exactly once. The factor 1/D\ accounts for the volume of a D-simplex. We now 
multiply the nominator and the denominator by e^("~i). Together with the factors ^/q^~^{v) 
and the factor from the integral, we get e^"'/V{v, e)"^"*^. The volumes in the denominator are 
absorbed into the Poisson brackets by the standard technique. The factors of Aaij are turned 
into holonomies {hs^)KL = Skl + ee"(0)Aa/j (t^'^)^^ + O(e^) using the the same amount of 
factors of e since we note that the zeroth order of the expansion of the holonomies vanishes when 
inserted into the Poisson brackets. We abbreviated = t, to simplify notation. 

The field strength tensors can be dealt with as follows. Let e, e' be arbitrary paths which 
are images of the interval [0, 1] under the corresponding embeddings, which we also denote by 
e, e' such that v = e(0) = e'(0). For any < e < 1 set e^{t) := e{et) for t G [0, 1] and likewise 
for e'. Then we expand he^{A) in powers of e. Consider the loop cxee,e'^ where in a coordinate 
neighbourhood 



ee(4t) 0<f<l/4 
e,{l) + eiiAt-l)-v l/4<t<l/2 
e'^{l) + e,{3 - 4t) - V 1/2 < t < 3/4 



e',(4-4t) 3/4<t<l. 

Now expanding again in powers of e we easily find ha , = lD+i + ^'^Fabi,jT^'^e"'{0)e"'{0) + O{e^). 
Since the indices of the field strength tensors are contracted only with other antisymmetric index 
pairs, the zeroth order of the expansion vanishes as well as the orders beyond in the limit e ^ 0. 
The remaining factors of e are absorbed into covariant derivatives using the approximation 



(/ie(0,e)7r«(e(e))/ie(0,e)-i-7r»(i;))''^ 



= ((1 + ee''(0)^)(7r^(^;) + ee%0)dcAvm ' ee''(0)^d) - Tr'iv)^ + 0{e^) 
= ee'^(0)r»c7r"^^(i;) + 0(e2). (5.17) 

We note that partial derivatives can be dealt with in the same way. 

At this point, all factors of e have been absorbed into holonomies and derivatives. It is 
key that the volume operators are ordered to the right in the quantum theory since then, the 
Hamiltonian constraint evaluated on a cylindrical function fy will only act on the vertices of 7. 
The action at vertices however does not depend on the value of e > and we can take the limit 
e ^ 0, thus removing the regulator. 

In order to quantise the Hamiltonian constraint, we have to replace the holonomies by 
multiplication operators, the volumes by volume operators, and the Poisson brackets hy i/h 
times the commutator. 



5.5 Reguleirised Quantities 

In order to construct a well defined Hamiltonian constraint operator, we have to express it in 
terms of operators well defined on the kinematical Hilbert space. Instead of writing down the 
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explicit regularisation for the proposed Hamiltonian constraint, we want to provide a toolkit 
for a general class of operators. In the following, we will propose "regulated" versions of the 
phase space variables, marked by an upper ^ in front. The idea will be to replace all phase space 
variables in the classical Hamiltonian constraint by their corresponding regulated versions, do 
some additional minor modifications and directly arrive at the Hamiltonian constraint operator, 
without explicitly dealing with the triangulation and the correct powers of e. Since the final 
constraint operator will only act on vertices of 7, it is sufficient to regularise the phase space 
variables at vertices 

In what follows, wc use a graph adapted coordinate system, meaning that the spatial coor- 
dinates a,b, . . . = 1, . . . , D enumerate the D edges incident at v of a D-simplex. 

5.5.1 L> + 1 > 3 Arbitrary 

We will express all the basic variables in terms of holonomies living on the edges of the adapted 
triangulation and volume operators acting on it. First, we notice that 

^(V^^+Va/j(^)) := ^-^^{hsJi''{{hsJ-K'j, {V{v, e)r+'} (5.18) 

is gauge covariant and reduces to ey/q^'^^'iTaij{v) in the limit e ^ 0. The factor of e is expected 
as the regulated quantity has a lower spatial index. In the end, when the complete constraint 
operator will be assembled, all factors of e will cancel out. We restrict x > —1 because powers 
of the volume operator will be defined by the spectral theorem in the quantum theory. 
For the KKEE terms, we propose 

X ^(^- Vl^^(^))(/i.jM^ { (/^ejpjv, {^i^WK e), e)}} , 

(5.19) 

where the '"tt"-^'^ will be defined below. 

Next, we regulate the gauge unfixing term DF~^D with density weight 1. We will place zero 
density into and a density weight of 1/2 into each D. Accordingly, 

/ 77-l\ iV M ^, rzA^EFGHN ^ I T?-^\ ^ ABCDM /'c on^ 

jcd,a6=7V^e TT(c\Ef{F )d)GH,{a\AB''b)CDe (5.20) 

becomes 

[VQ ^ )cd,ab=l(- W1T^{c\Ef) [y/q t )a)GH{a\AB ^ ' 

with 

^(VgV-i)^^^^^^^ := ^^J_-^^(^^„^^)^(^v^,BD)(^(^/^-V^^^^y(V^v^ei,^)-7?^^) 

(^^AB^K[I^J]L _ 2^^^^S[/^J]X^ _ (5 22) 

The D constraint contains a covariant derivative which we regularise as 

^{^-'DaTT'^^) := (hsJ{VQ-'Asa))h7: ' ^(^^^ V'(^))) • (5-23) 
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The full D constraint 



= -euKLM-"'' (-('^I^^I^c^^) V) (5.24) 

can thus be regularised as 

^^-'^'D'^) = -ejjKLj^'iV^-'^^^'') (^(V^-V('^l^^)^(Vg"'l)cVr^)V)) . (5.25) 

A different regularisation procedure for the DF^^D part of the Hamiltonian constraint which 
is based on field strength tensors is outlined in appendix [A} 

In general, a generic power of l/^/q needed to turn the individual terms with densities > 1 
into densities of weight 1 can be constructed as 

J 1 \ /l^^ 



{-2xD-2) / {2 



detl^i^^^^'TTaun^'^W')) (5.26) 



with the usual x > —1. 

The field strength tensors are regularised as 



'FaMJ={ha.^..^)^jl^S^^ (5.27) 

while we set 

^{A^iAv), ■} = -(/i.J/^{(/i.r'kj, •}• (5.28) 

5.5.2 D + 1 Even 

Let n = {D- l)/2. We "regulate" 

'{^/q^^~^'^''7r''^■^{v)) ^ ^ e"^^^^-''"""e^'^-^i-^^-^-^"sgn(dete)(^) 



(D-l) 

(V^('+^V,,,,^,(^;))^(Vg(i+^V,,^,^H^)) • • ■ 
{V<l^'-''^^b„i.KAv)nVQ^'''^^^c„J.''"iv)) (5.29) 



and 



^(Sgn(det(ei))) « ^,IJhJ^...InJ.^aa,b,...a„b„e^^{D~l)/D^^^^^ 

^(v^^^-^)/^^....i.j^(V^(^-^)/^n,..^^)... 

^(V^(^-^)/^vr.„,„^J^(V^(^-^)/^.,„,„^"). (5.30) 
For the Euclidean part of the Hamiltonian constraint, we need 

/ [a\IK b]J \ 1 
e I ^ \ ^ ^ abcaibi...an~ib„-i IJKLhJi...I„-iJ„-i (J \ 

^{VqT^cKLyi^/q-^aihKiYi^/q'^biJi^'^) ■ ■ ■ 

^(^/^vra„_i/„_ii^„_J^(V^vrfe„.,j_,^"-). (5.31) 

As stressed before, the two possibilities to express the Euclidean part of the Hamiltonian con- 
straint are equally complicated. 
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5.5.3 D + l Odd 

Let n = {D- 2)/2. We "regulate" 

^(V^(^-i)^'7r'^^-^(^;)) « ^^^^^, e"''''i^^---^"^"e^^^^^^^---^"-^"sgn(dete)(^)^(Vg^^+"VfeLi^(i;))^n^(^) 



and 



For the Euclidean part of the Hamiltonian constraint, we need 



'[a\IK^b].j\ 1 



\ ^ J 2{D-2)\ 
5.6 The Hamiltonian Constraint Operator 

At this point, we are ready to assemble the Hamiltonian constraint operator. The general idea 
of the regularisation has been described in section 5.4 Here, we provide a toolkit in order to 
assemble the constraint operator. 

(1) The Euclidean part %e = 2^''^"^^^^^'^ ^Fabij of the Hamiltonian constraint can be quan- 
tised with the methods described above and using the following recipe. The corresponding 
operator can then be used in commutators to express additional parts of the full Hamilto- 
nian constraint operator. 

(2) Use classical identities in order to express the Hamiltonian constraint in terms of connec- 
tions Aaij, volumes V{x,e) and Euclidean Hamiltonian constraints '?^£;(x,e). 

(3) Replace all phase space variables by their corresponding regulated quantities. 

(4) Instead of the the integration d^x, put a sum J2v(^v{j) '^^^^ vertices v of the 
graph 7. 

(5) For every spatial e-symbol, put a sum -J^ Ylv{A)=v ^^^^ D-simplices having v as 
a vertex. The holonomies associated with the e-symbol are evaluated along the edges 
spanning A. 

(6) Substitute the Poisson brackets by ^ times the commutator of the corresponding operators, 
i.e. the multiplication operator he and the volume operator V. 

In order to understand the double sum over D-simplices appearing in the KKEE and the 
gauge unfixing term, consider the following argument given in a similar form in |22j : Since 
limg^o {^/^^)Xe{x,y) = S^{x,y) we have limg^o ^) = y/(l{x) ■ It is also easy to see that 
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' {A,KL{y),v{y,e)}z'^''Hy) 



for each e > we have that dV/dir"'^'^ (x) = 6V{x,e)/6'iT"'^'^ (x). The terms under consideration 
are of the form 

,D Vqix)7Tau{x)Z-''ix)^{x)7ThKL{x)Z'^^{x) 

d x^ , 5.35 

where Z""^^ is a density of weight +1 and stands symbolicahy for the remaining terms, including 
a spatial e-symbol with upper indices, one of which is a. We rewrite this expression as 

;^o^'(^-^) I ' " — — J ^ — — 

= lim 4(25- 1)2 y d^x{A,,j(x),yy(^;^}Z«^-^(x) y" d^yXe(x,y){A,,^i.(y),yy(^}Z^^^(y). 

Triangulation leads to two sums over vertices and two sums over D-simplices containing the 
individual vertices. In the limit e — )■ however the two sums over vertices collapse to a single 
sum over vertices due to the Xe term and we have the desired result. 

5.7 Solution of the Hamiltonian Constraint 

As in the 3 + 1-dimensional treatment, we realise that the only spin changing operation of 
the Hamiltonian constraint is performed by its Euclidean part. The construction of a set of 
rigorously defined solutions to the diffeomorphism and the Hamiltonian constraint described in 
|41j thus immediately generalises to our case. 

5.8 Master Constraint 

The implementation of the Master constraint 

M=-/'d^x^%4 (5-37) 
2y. ^{x) ^ ' 

works analogously to the 3 + 1-dimensional case described in ||42j. The inverse square root is 
split up between the two Hamiltonian constraints and hidden by adjusting the power of the 
volume operators as before. The result of the derivation is the Master constraint operator 

MTy,^ :=^Qm(T[,,],T[,])T[,,] (5.38) 

with 



Qm{1,1') = Y.i1[s] J2 liClTs,iis]))l'{ClT,^as])) (5.39) 

[s] ■oGy(7(soW)) 

and /(C'lr^Q([<,])) being the evaluation of / on the Hamiltonian constraint operator with the 
additional 1 / ^ hidden in the volume operator (s) . The proof of the following theorem generalises 
with obvious modifications from the treatment in |3]. 
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Theorem. 



(i) The positive quadratic form Qm is closable and induces a unique, positive self-adjoint 
operator M on Hdiff- 

(a) Moreover, the point zero is contained in the point spectrum of M . 

We deal with the problem of Tidis not being separable by using ^-equivalence classes of 
spin-networks, see [32] . Now, a direct integral decomposition of T-L^^^ is available: 

Theorem. 

There is a unitary operator V such that VT-L^^^g is the direct integral Hilbert space 

r df,{x)n%ff{x) (5.40) 

where the measure class of jjL and the Hilbert space ^^^■^(A), in which VMV~^ acts by multipli- 
cation by X, are uniquely determined. 

The physical Hilbert space is given by H^^^yg = ?^^j^(0) . 

We notice that we could define an extended Master Constraint that also involves the sim- 
plicity constraint. 

5.9 Factor Ordering 

In |34| I35j . it has been shown that there is a unique factor ordering which results in a non- 
vanishing flux operator expressed through the volume operator and holonomies in the usual 
3 + 1 dimensional LQG. The idea, translated to our case, is that the volume operator in the 
expression for '^■k"-^'^ has to act on an at least L'-valent non-planar vertex and the holonomies 
in the expression have to be ordered to the right for this to be ensured. Apart from ordering 
individual terms of the sums appearing differently (which would be highly unnatural), this leaves 
only one possible factor ordering. We remark that the proof of the equivalence of the "normal" 
and "derived" flux operator given in |341 [35] does not generalise trivially to our case since it 
is explicitly based on SU(2) as the internal gauge group. We leave this point open for further 
research. 

In order to ensure that the Hamiltonian constraint only acts on vertices, we order in all three 
terms either a commutator [/i^^,t^] or a double-commutator [h~^^ '(He,V] to the right. 

We leave the remaining details of the factor ordering open, as this paper only intends to 
show that a quantisation is possible in principle. 

5.10 Outlook on Consistency Checks 

At this point, one might ask if there are good indications whether the proposed theory is phys- 
ically viable. In case of the usual formulation of LQG in terms of Ashtekar-Barbero variables, 
it was shown in [43j that a quantisation of Euclidean General Relativity in three dimensions 
with methods very similar to the ones used in LQG recovers the known solutions of three- 
dimensional General Relativity familiar from other approaches. The reason why these theories 
match is that they both use the gauge groups SU(2) and that a suitable redefinition of the 
Lagrange multipliers of Euclidean three-dimensional General Relativity leads to a Hamiltonian 
constraint with the same algebraic structure as the Euclidean part of the constraint familiar from 
LQG. A similar check is conceivable for the presented theory in that we can describe Lorentzian 
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three-dimensional General Relativity using SU(2) as a gauge group, which would result in a 
different Hamiltonian constraint. One could now check if the solution space of Lorentzian three- 
dimensional General Relativity is reproduced when using SU(2) as a gauge group and thus 
mimicking the internal signature switch which is also done in this formulation. 

As for the simplicity constraint, we cannot use three-dimensional General Relativity as a 
testbed since the simplicity constraints only appear in four and higher dimensions. In this 
paper, two different regularisations for the gauge unfixing part of the Hamiltonian constraint 
were introduced, one in section 5.5 and one in appendix|Xj While the regularisation introduced 
in section 5.5 preserves the closure of the quantum constraint algebra, this is not obvious for 
the regularisation in appendix [A] since terms quadratic in the field strength appear. 

Another approach to consistency checks is to compare our formulation in four dimensions 
to the usual LQG formulation. In section 4.1, the area operator was shown to have the same 
spectrum as in standard LQG, which however does not come as a surprise regarding similar 
results from spin foam models. As for the volume operator, we do not know whether the spectrum 
matches the one of standard LQG. This is also tied to the fact that we are only interested in 
the spectrum on the solution space to the vertex simplicity constraint operators, for which we 
do not have a completely satisfactory proposal. We remark that a matching spectrum of the 
volume operator can be obtained by using a weak implementation of the linear vertex simplicity 
constraints [Hj. However, as explained in our companion paper [25], this approach comes with 
its own problems in the canonical theory. 



6 Conclusion 

In this paper we have demonstrated that by a straightforward adaption of the toolbox developed 
for LQG in 3 -|- 1 dimensions also the constraints of our new connection formulation of General 
Relativity in any dimension Z? -|- 1 > 3 can be quantised analogously and rigorously. The higher 
dimension does not require much more complexity than in 3 -|- 1 dimensions. We conclude 
that our new connection formulation has a consistent quantisation. The next task is to study 
matter coupling, in particular coupling to supersymmetric matter in interesting dimensions, 
where String theories and Supergravity theories are defined, and the quantisation thereof. This 
has to be done, as in 3 -|- 1 dimensions, in a background independent way, a task to which we 
turn in the next papers of this series [45 1 H6l HT] . 

In four dimensions, we now have the special situation that there are two formulations of 
LQG, one based on the usual Asthekar-Barbero variables, and one based on the variables pro- 
posed in this series of papers. From a direct comparison, one concludes that the new formulation 
is more complicated since the Hamiltonian constraint contains an additional term resulting from 
gauge unfixing. Two different regularisations for this term were introduced, the first one directly 
regularises the covariant derivatives in this term, the second one uses a Poisson bracket identity 
involving the Field strength and the whole expression is thus quadratic in the field strength. 
Both of these regularisations do no appear in the standard case and the Hamiltonian constraint 
operator is thus more complicated. On the other hand, since it is already hard to deal with 
the usual Hamiltonian constraint, we cannot conclude that our Hamiltonian constraint is sig- 
nificantly more complicated. The main problem remains the simplicity constraint for which a 
satisfactory implementation has to be found which is compatible with the action of the Hamil- 
tonian constraint and allows for a unitary map to the Ashtekar-Lewandowski Hilbert space. 
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A Alternative Regular isat ion of DF 

It was suggested by Wieland [H] that one could simplify the D constraints by using the classical 
identity 

2D[,^7T,]ij{x) = -{D - l){Fabij{x), V{x, e)}, (A.l) 

i.e. the torsion of the gravitational connection can be expressed using a Poisson bracket which 
will become a commutator in the quantum theory. Since the D constraints appear quadrati- 
cally in Ti, this type of regularisation results in a more non-local operation of the Hamiltonian 
constraint. 

In order to apply the above identity, we recall from [2j that we can extend the covariant 
derivative in 

= -^IJKLM-''' {.(^^^^D,."^^^) (A.2) 
by a Christoffel symbol acting on spatial indices on the constraint surface. Therefore, 

K = -^IJKLM^''' (vr('^l^^gg'')^Z),vr/^) (A.3) 

and we calculate 



-(i-{a\ABT^h)CD^ ^IJKLM^ 11 ^[c^d] N ^ 

^ {D- 3)!(Z) - ^'^^dbKLiF'f^-^'^^^ (A.4) 

with 

^pyU,bKL ^ 2 (^-E-\K]Eb[I^J][L + qq-b^^HflL^ . (A.5) 
In order to have direct access to K^Jjj^ we can invert F' as 

{F' ^)bKL,cMN = (^-^qbcVMlKVL]N + 2^b\M]Ec[KVL]lN^ (A. 6) 
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and write 



^ ^ ^ 8((D-3)!)2(D-1)2 ^^-^^ 



^mAByj" )a)CD,eMN^ ) fOP,{c\EF^d)GH(^ ]V 



8((D-3)!)2(D-1)2 



with 



j^/abABCDA-f (^r, ^5 -l'^P \ .jr ^EFGHN j^icd 



^M = -^IJKLM^ m'U[c'^d]N]- (A 



We can also implement the above Poisson bracket identity without starting from the original 
D constraints but by trying to find an easier expression for % —% directly from D[j,7r^/j. It 
turns out that 

■H-H^ aD[,VQ^,^ljm,VQ^,]KLWn'^ [<i"'E-^E'' + \qq''^' q'^'^W^ ■ (A.9) 

The obvious question at this point is which of the two expressions is suited better for a 
quantisation. Although a satisfactory answer might only be possible after studying the quan- 
tum dynamics, we see at the classical level that the second expression has a less complicated 
index structure due to the missing epsilon symbols. On the other hand, it contains correction 
terms proportional to K^j , which are absent due to the epsilon symbols in the first expression. 
In the formulation studied in this paper, this does not affect the theory since K'f ~ on the 
constraint surface In general, this won't be true any more when coupling fermions |45] or 
using the time normal as a independent field |46j in other papers of this series. Although 
introducing additional correction terms, an independent time normal would simplify the ex- 
pression since the action of a multiplication operator corresponding to is simpler than the 
regularised version of n^n^ [i:). 
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